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Abstract 

The  dynamic  stall  flow  field  over  a  NACA  0012 
airfoil  pitching  transiently  from  0  —  60®  at  a  con¬ 
stant  rate  under  compressible  flow  conditions  has 
been  studied  using  the  real-time  technique  of  point 
diffraction  interferometry.  This  investigation  using 
nonintrusive  diagnostics  provides  a  quantitative  de¬ 
scription  of  the  overall  flow  field,  including  the  finer 
details  of  dynamic  stall  vortex  formation,  growth  and 
the  concomitant  changes  in  the  pressure  distribution. 
Analysis  of  several  hundred  interferograms  obtained 
for  a  range  of  flow  conditions  shows  that  the  peak 
leading  edge  suction  pressure  coefficient  at  stall  is 
nearly  constant  for  a  given  free  stream  Mach  number 
at  all  nondimensional  pitch  rates.  Also,  this  value  is 
below  that  seen  in  steady  flow  at  static  stall  for  the 
same  Mach  number,  indicating  that  dynamic  effects 
significantly  affect  the  separation  behavior.  Further, 
for  a  given  Mach  number,  the  dynamic  stall  vortex 
seems  to  form  rapidly  at  nearly  the  same  angle  of 
attack  for  all  pitch  rates  studied.  As  the  vortex  is 
shed,  it  induces  an  anti-clockwise  trailing  edge  vor¬ 
tex,  which  grows  in  a  manner  similar  to  that  of  a 
starting  vortex.  The  measured  peak  suction  pressure 
coefficient  drops  as  the  free  stream  Mach  number  in¬ 
creases.  For  free  stream  Mach  numbers  above  0.4, 
small  multiple  shocks  appear  near  the  leading  edge 
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above  the  shear  layer.  The  unsteady  flow  peak  suction 
pressure  coefficient  at  any  angle  has  been  experimen¬ 
tally  demonstrated  to  lag  that  of  the  steady  flow  at 
the  same  angle  of  attack. 

Nomenclature 

Cp  pressure  coefficient 

c  airfoil  chord 

M  free  stream  Mach  number 

Re  Reynolds  number  based  on  c  and  Uoo 

Uoo  free  stream  velocity 

x,y  chordwise  and  vertical  distance 

a  angle  of  attack 

a  pitch  rate  in  deg/sec 

a"''  nondimensional  pitch  rate, 

7  ratio  of  specific  heats 

e  fringe  number 

p  density 

pr  density  at  reference(atmospheric)  conditions 

1.  Introduction 

The  utilization  of  dynamic  stall  as  a  method  for 
increasing  the  maneuverability  and  agility  of  aircraft 
has  received  significant  attention  during  the  past  few 
years.  Several  researchers*"®  have  studied  the  flow 
over  pitching  airfoils  using  flow  visualization  and  un¬ 
steady  pressure  measurements  and  provided  valuable 
information  on  the  dynamic  stall  phenomenon.  How¬ 
ever,  all  these  studies  are  at  low  speeds.  Lorber  and 
Carta'  have  obtained  measurements  under  compress¬ 
ibility  conditions  which  showed  that  at  the  higher 
Mach  number,  the  flow  could  not  develop  the  suc¬ 
tion  levels  observed  under  incompressible  conditions, 
and  noted  that  this  effectively  limited  the  stall  delay 
that  could  be  achieved.  The  computational  studies  of 
Visbal*  lend  some  support  to  this  result.  Since  com¬ 
pressibility  effects  have  been  shown®  *®  to  change  the 
way  that  dynami'-  stall  develops,  a  better  understand¬ 
ing  of  these  effects  hcis  been  of  interest  in  the  develop¬ 
ment  of  supermaneuverable  aircraft  and  highly  agile 
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encoder,  whose  output  was  input  to  the  digital  I/O 
board  of  a  micro  VAX  II  Work  Station  and  timed  with 
its  internal  clock.  Figure  2  shows  an  example  of  the 
actual  rates  obtained,  including  the  variation  of  the 
angle  of  attack  during  the  various  parts  of  the  pitch- 
up  motion  in  an  experiment.  As  can  be  seen,  the 
airfoil  goes  through  its  static  stall  angle  with  a  linear 
rate  of  change  of  angle  of  attack.  For  the  highest  rate, 
the  motion  is  completed  in  18  ms,  beyond  which  the 
system  is  settling  down  (at  the  highest  angle).  All  the 
tests  were  limited  to  the  linear  range. 

The  CDSF  is  equipped  with  a  wide  range  of  non- 
intrusive  optical  flow  diagnostic  instrumentation  such 
as  stroboscopic  schlieren,  laser  Doppler  velocimetry, 
holographic  interferometry  and  point  diffraction  in¬ 
terferometry  systems.  The  present  paper  will  present 
results  obtained  using  the  point  diffraction  interfer¬ 
ometry  system. 

2.3.  The  Point  Diffraction  Interferometry 
Technique 

The  point  diffraction  interferometry  technique 
used  in  this  study  utilizes  the  ability  of  a  point  dis¬ 
continuity  (in  the  form  of  a  pin-hole)  located  at  the 
image  of  a  point  source  to  diffract  a  portion  of  the 
incident  light  into  a  spherical  reference  wave  front. 
In  the  present  application,  the  primary  optics  of  an 
existing  schlieren  system  were  used  (see  Reference  13 
for  details),  with  a  pulsed  Nd:YAG  laser  replacing 
the  conventional  spark  as  the  light  source,  and  a  spe¬ 
cially  created  point  diffractor  replacing  the  usual  knife 
.edge.  The  laser  light  was  expanded  through  a  micro¬ 
scope  objective  to  fill  the  schlieren  mirror,  transmit¬ 
ted  through  the  test  section,  and  refocused  by  another 
schlieren  mirror.  An  exposed  photographic  plate  used 
to  create  the  point-diffraction  spot  was  placed  at  the 
focus  of  this  second  mirror,  and  the  laser  was  pulsed 
with  enough  energy  to  burn  a  hole,  or  spot,  in  the 
emulsion  located  at  the  focal  plane  of  the  second  mir¬ 
ror.  The  spot  was  created  in  sitti  by  passing  light 
through  the  test  section  at  a  no-flow  condition.  The 
spot  was  precisely  tailored  to  the  application  under 
investigation,  automatically  correcting  for  nonunifor¬ 
mities  in  the  light  source  or  optics.  With  the  flow 
turned  on,  the  laser  was  triggered  externally  at  the 
desired  angles  of  attack  and  the  real-time  interference 
fringes  were  recorded  on  Polaroid  film(ASA  3000), 
which  were  available  for  immediate  viewing.  More 
details  about  the  PDI  technique  can  be  found  in  Ref¬ 
erences  13  and  14. 

2.4.  Interferogram  Image  Processing 

The  quantitative  nature  of  the  interferograms  en¬ 
abled  computation  of  the  pressure  distribution  over 
the  airfoil  when  the  flow  was  attached.  The  inter¬ 
ferograms  were  processed  in  a  semi-automatic  mode 
using  a  specially  developed  software  package.  The 
program  read  a  digitized  (256  gray  level)  interfero¬ 
gram  on  an  IRIS  Work  Station  and  overlayed  an  air¬ 
foil  using  the  triangular  registration  markers  seen  in 
the  images.  The  intersection  of  the  fringes  with  the 
airfoil  upper  a»;d  lower  surfaces  (or  the  local  boundary 
layer  edge,  when  detectable)  were  interactively  picked 
by  the  user.  Since  each  fringe  is  a  line  of  constant  den¬ 
sity,  the  corresponding  pressure  at  the  boundary  layer 


edge  could  be  derived  using  isentropic  flow  relations. 
This  pressure  was  then  used  as  the  surface  pressure, 
under  the  boundary  layer  assumptions.  The  density 
along  any  fringe  could  be  calculated  from  the  Glad¬ 
stone  -  Dale  equation*^,  which  for  the  present  wind 
tunnel  and  laser  simplifies  to 

p  —  pr  =  0.00942l€ 

where  e  the  fringe  number  is  0,  ±1,±2,...  for  the 
bright  fringes  and  ±5,  ±|,  ±|,...  for  the  dark  fringes. 
Fringes  from  the  free  stream  to  the  stagnation  point 
have  positive  values.  The  corresponding  Cp  distribu¬ 
tions  were  then  computed  from  the  relation 


In  the  end,  the  package  provided  an  output  data  file 
containing  the  various  physical  variables,  in  a  format 
suitable  for  plotting.  Typical  processing  time  was 
about  3-5  minutes  per  image. 

In  cases  where  the  fringe  density  was  high  or  the 
fringes  were  fuzzy,  the  user  could  go  into  the  ‘off- body’ 
mode  and  pick  fringes  along  a  line  parallel  to  and 
away  from  the  airfoil  surface  where  the  fringes  are  far¬ 
ther  apart.  For  this  purpose,  an  option  to  superpose 
two  larger  airfoils  over  the  image  on  the  screen  was 
provided.  The  fringe  intersections  on  the  larger  air¬ 
foils  were  then  projected  to  the  airfoil  surface  suitably. 
At  angles  of  attack  near  the  dynamic  stall  angle,  the 
fringes  near  the  leading  edge  region  were  very  dense 
reflecting  the  large  local  density  gradients.  Further, 
in  this  region,  optical  noise  introduced  by  the  shad¬ 
owgraph  effect  generally  lowered  the  contrast,  making 
it  a  location  where  the  off-body  mode  needed  to  be 
invoked. 

In  the  present  study  the  entropy  change  in  the 
vortical  flow  was  ignored  (for  lack  of  a  better  method). 
Interferograms  with  shocks  have  not  been  processed 
because  of  this  limitation. 

2.5.  Experimental  Conditions 

The  Mach  number  of  the  experiment  ranged  from 
0.2  -  0.45.  The  corresponding  nondimensional  pitch 
rate  varied  from  0.020  to  0.045,  depending  on  the 
Mach  number.  The  actual  maximum  rate,  as  stated, 
was  limited  to  3600  degrees/sec.  The  images  were 
obtained  at  a  resolution  of  0.5  degrees,  with  some  ex¬ 
ceptions  when  the  results  indicated  need  for  a  better 
resolution. 

The  nondimensional  pitch  rates  reported  are 
based  on  the  total  time  for  pitching  from  0°  —  57°. 
However,  the  hydraulic  control  system  of  the  pitching 
drive  caused  the  airfoil  to  pitch  5  -  8%  faster  in  the 
0°  —  10°  range  and  sometimes  in  the  0°  —  30°  range. 
This  difference  is  not  believed  to  significantly  affect 
the  global  results  of  the  study. 

The  experiments  were  conducted  in  two  phases. 
In  the  first  phase,  the  full  flow  field  interferograms 
were  obtained  for  a  range  of  conditions;  phase  II  fo¬ 
cused  on  the  leading  edge  flow  details  only. 
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3.  Results  and  Discussion 

The  interferogram  images  will  be  first  discussed 
qualitatively  as  flow  visualization  images.  The  quan¬ 
titative  pressure  distributions  derived  from  the  images 
will  be  presented  in  the  second  part  of  the  section. 
Whereas  much  of  the  paper  is  devoted  to  the  leading 
edge  flow,  when  appropriate,  the  full  flow  field  is  also 
discussed. 

3.1.  Discussion  of  Interferogram  Images 

A.  Separation  Bubble  and  Dynamic  Stall 

Fig.  3  is  a  point  diffraction  interferogram  of  the 
flow  at  M  =  0.3,  a  =  12°  and  a'*’  =  0.03.  This  im¬ 
age  reveals  some  important  features  of  the  flow.  The 
dark  closed  fringe  on  the  lower  surface  slightly  aft  of 
the  leading  edge  is  the  stagnation  point  fringe.  The 
suction  pressure  developed  by  the  airfoil  causes  the 
local  flow  to  accelerate,  resulting  in  strong  density 
changes,  which  is  seen  in  the  figure  as  a  concentration 
of  fringes  near  the  leading  edge  on  the  upper  surface. 
The  close  spacing  of  the  fringes  also  means  that  the 
flow  gradients  are  very  high.  In  fact,  21  dark  fringes 
are  present  in  about  1  millimeter  (x/c  sa  0.1)  in  this 
image,  indicating  that  the  local  maximum  Mach  num¬ 
ber  is  0.71  and  the  local  pressure  coefficient  is -3.75  at 
x/c  =  0.01.  Downstream  of  this  point,  a  steep  adverse 
pressure  gradient  region  develops(see  Sec.  3.2.C)  and 
classical  laminar  flow  separation  occurs.  The  sepa¬ 
rated  shear  layer  reattaches  after  it  transitions  into  a 
turbulent  layer,  forming  a  bubble.  In  the  figure,  the 
fringes  in  the  bubble  appear  as  lines  emanating  from 
the  leading  edge  which  turn  abruptly  towards  the  up¬ 
per  surface  at  x/c  fs  0.05  and  turn  sharply  again  as 
the  local  boundary  layer  is  encountered.  Inside  the 
bubble,  the  fringes  run  parallel  to  the  surface  locally, 
representing  the  pressure  plateau  normally  associated 
with  laminar  separation  bubbles.  The  accompany¬ 
ing  pressure  distribution  (which  will  be  discussed  in 
Sec.  3.2. C)  obtained  by  processing  the  image  using 
the  fringe  analysis  software  shows  the  suction  peak, 
the  drop  in  suction  due  to  the  adverse  pressure  fol¬ 
lowing  it  and  the  laminar  separation  bubble,  which 
is  indicated  by  the  plateau  in  the  distribution.  In 
the  interferograms  to  be  discussed,  all  of  the  features 
upstream  of  the  bubble  are  found.  As  dynamic  stall 
occurs,  differences  evolve  which  will  be  pointed  out. 

B.  Leading  Edge  Supersonic  Flow 

Chandrasekhara  et  al*®  found  that  the  airfoil 
leading  edge  flow  can  become  supersonic  leading  to 
formation  of  multiple  shocks.  Fig.  4  provides  quan¬ 
titative  documentation  of  this.  Depending  upon  the 
angle  of  attack,  a  shock,  or  multiple  shocks  form  in 
the  flow.  Fig.  4  shows  a  PDI  image  for  M  =  0.45, 
at  o  =  12.6°  and  or"^  =  0.0313.  Fringe  counting 
shows  that  the  local  Mach  number  abpad  of  the  shock 
is  greater  than  1.0  and  at  the  foot  of  the  shock,  it 
is  about  1.2.  Although  the  flow  is  only  weakly  su¬ 
personic,  the  shock  causes  the  leading  edge  laminar 
boundary  layer  to  separate*®.  This  separated  free 
shear  layer  develops  a  waviress,  which  causes  the  flow 
downstream  of  the  shock  to  go  through  a  series  of  ac¬ 


celerations  and  decelerations.  As  the  flow  negotiates 
the  crests  and  valleys  of  this  wavy  shear  layer,  expan¬ 
sion  waves  and  compression  waves  develop,  causing  a 
series  of  shocks.  The  latst  shock  in  the  series  appears 
to  be  the  strongest  and  the  flow  becomes  subsonic 
downstream.  The  occurrence  of  multiple  shocks  is  re¬ 
peatable,  and  the  shocks  were  found  to  be  present 
over  an  angle  of  attack  range  of  about  one  degree. 

C.  Flow  Description  at  M  =  0.2,  o'*"  =  0.03 

Fig.  5  presents  a  montage  of  interferograms  for 
M  =  0.2  at  a  nondimensional  pitch  rate  of  0.03  for 
12°  <  a  <  24.75°.  A  separation  bubble  is  present 
for  Q  =  12°;  analysis  of  the  interferograms  showed 
that  the  bubble  first  appeared  at  a  =  7°  (as  opposed 
to  about  6  degrees  in  steady  flow  at  M  =  0.2).  At 
o  =  15°,  the  fringe  near  the  trailing  edge  indicates  a 
mild  local  flow  separation  as  this  fringe  first  moves 
into  the  wake  and  turns  sharply  back  towards  the 
trailing  edge.  Also,  the  leading  edge  bubble  starts 
to  open-up  and  the  first  imprint  of  the  dynamic  stall 
vortex  becomes  distinct  at  a  =  17°.  As  it  grows,  the 
enveloping  shear  layer  moves  downstream  and  the  air¬ 
foil  boundary  layer  thickens.  The  growth  of  the  vor¬ 
tex  continues  while  it  convects  over  the  airfoil  dur¬ 
ing  the  continuous  ramping  motion.  For  example,  at 
a  =  18°,  the  center  of  the  vortex  is  at  x/c  =  0.1, 
but  the  downstream  edge  of  the  surrounding  shear 
layer  is  at  x/c  =  0.3.  The  number  of  fringes  is  seen 
to  increase  till  a  =  18°  demonstrating  that  the  low 
pressure  region  continues  to  grow.  By  a  =  20°,  the 
leading  edge  flow  has  separated  and  the  fluid  aft  of 
the  shear  layer  enclosing  the  vortex  has  been  con- 
vected  past  the  trailing  ^ge  and  thus,  no  fringes  can 
be  seen  in  this  region.  Once  the  vortex  grows  and 
begins  to  convect,  the  innermost  fringes  become  cir¬ 
cuit  (for  example  at  a  =  21°  and  22°),  whereas  the 
outer  fringes  still  end  on  the  surface.  The  number  of 
fringes  inside  the  vortex  is  now  larger  than  that  at 
the  leading  edge.  This  implies  that  the  pressure  in 
the  vortex  core  is  lower  than  the  peak  suction  over 
the  airfoil.  Surface  pressure  measurements,  however, 
cannot  reveal  this.  The  vortex  is  eventually  shed  by 
a  =  24.75°,  and  the  flow  reaches  the  deep  stall  state. 
A  counter  clockwise  trmling  edge  vortex  also  forms  at 
this  angle  of  attack.  For  this  test  condition,  the  vor¬ 
tex  remains  over  the  surface  for  a  large  angle  of  attack 
range  of  12  degrees.  Similar  studies  on  an  oscillating 
airfoil  (o  =  10°  -1-  10°sin  wt)  at  M  =  0.2  showed  that 
deep  stall  occurred  at  q  =  18.1°.  The  presence  of  the 
dynamic  stall  vortex  till  a  =  24.75°  elicits  the  fact 
that  motion  history  plays  a  key  role  in  the  dynamic 
stall  process.  In  this  case,  the  transiently  pitching 
airfoil  is  found  to  be  better  in  sustaining  the  dynamic 
lift  generated. 

D.  Flow  Sequence  at  M  =  0.45,  =  0.02. 

Fig.  6  shows  a  similar  set  of  interferograms  for  M 
=  0.45  at  a  pitch  rate  of  0.02.  At  thio  Mach  uumber, 
compressibility  effects  dominate*®.  The  large  number 
of  fringes  seen  at  low  angles  of  attack  are  due  to  the 
larger  density  changes  in  the  flow  at  this  higher  Mach 
number.  The  innermost  closed  fringe  intersecting  the 
lower  surface  near  the  leading  edge  encloses  the  stag¬ 
nation  point.  There  are  58  fringes  (corresponding  to 
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a  local  mach  number  of  ss  1.2)  around  the  leading 
edge  in  the  first  1.5%  of  the  airfoil  chord  at  o  =  9®. 
A  separation  bubble  forms  in  this  case  also  and  is 
seen  clearly  at  q  =  7®.  In  steady  flow  at  this  Mach 
number,  the  bubble  was  first  found  for  0  =  5®,  but  a 
delay  was  observed  in  its  first  appearance  in  the  un¬ 
steady  case.  At  Q  =  9®,  the  first  signs  of  the  dynamic 
stall  vortex/vortical  structure  are  seen  as  a  thin  shear 
layer  between  the  bubble  and  the  airfoil  upper  surface 
near  the  leading  edge  region(see  also  Fig.  7).  Mul¬ 
tiple  shocks  can  be  seen  above  the  airfoil  shear  layer 
and  the  flow  is  similar  to  that  discussed  in  Fig.  4. 
The  shocks  remain  on  the  surface  till  o  =  11®,  even 
after  the  dynamic  stall  process  is  well  under  way.  It  is 
not  yet  clear  whether  the  first  shock  induced  bound¬ 
ary  layer  separation  caused  the  dynamic  stall  vortex 
to  form. 

The  vortex  grows  as  the  shear  layer  enveloping  it 
reattaches  further  down  the  airfoil  towards  the  trail¬ 
ing  edge  with  increase  in  angle  of  attack.  As  in  the 
low  Mach  number  case,  the  fluid  aft  of  this  point  has 
been  swept  away  into  the  wake.  This  is  a  region  of 
nearly  stagnant  fluid.  No  circular  fringes  are  found 
in  the  vortex;  instead,  only  half  circular  fringes  are 
seen.  At  a  =  14.5®,  the  edge  of  the  shear  layer  has 
reached  the  trailing  edge,  and  deep  stall  occurs  at 
a  =  15.5®.  The  entire  sequence  lasts  only  6  -7  de¬ 
grees  in  angle  of  attack,  which  is  considerably  smaller 
than  the  12®  range  seen  for  M  =  0.2.  Although  the 
lower  Mach  number  result  discussed  was  obtained  at 
a  nondimensional  pitch  rate  which  was  50%  higher, 
the  range  of  angle  of  attack  over  which  dynamic  lift  is 
sustained  is  nearly  twice  that  seen  at  M  =  0.45.  This 
confirms  the  result  obtained  from  the  earlier  schlieren 
studies^®  that  compressibility  promotes  stall.  These 
results  agree  with  those  reported  in  Ref.  17  for  the  os¬ 
cillating  airfoil.  However,  it  should  be  noted  that  the 
two  different  motion  histories  will  force  differences  in 
the  details  of  separation.  Dynamic  stall  of  oscillating 
airfoils  is  influenced  by  hysteresis  which  affects  all  as¬ 
pects  of  the  flow.  The  transiently  pitching  airfoil  flow 
is  free  of  this  parameter.  In  addition,  since  the  airfoil 
continues  to  pitch  to  60®,  which  is  well  beyond  the 
static  stall  angle,  the  development  of  the  post-stall 
flow  and  the  interactions  at  the  trailing  edge  can  be 
studied  (discussed  in  Sec.  F),  which  is  of  importance 
to  the  supermaneuverability  problem.  This  informa¬ 
tion  is  also  of  value  in  comparing  computed  results 
against  experiments. 

The  differences  in  the  vortex  size  and  structure 
also  imply  that  the  overall  (global)  pressure  field  in 
the  incompressible  and  compressible  flow  cases  are 
quite  different.  The  sustained  presence  of  a  low  pres¬ 
sure  region  over  the  airfoil  upper  surface  is  evidence 
that  even  in  the  compressible  case,  there  is  increased 
lift  generated  during  the  dynamic  stall  process. 

E.  Leading  Edge  Flow  Stall  Vortex  Formation 

Carr  et  al  found  that  for  an  oscillating  airfoil, 
the  dynamic  stall  vortex  formed  just  as  the  separation 
bubble  burst.  It  has  now  been  found  to  be  true  for  the 
transiently  pitching  airfoil  also.  Fig.  7  presents  some 
interferograms  that  demonstrate  this  result  for  M  = 
0.3  and  q"*"  =  0.03.  At  a  =  15®,  the  fringes  enclos¬ 
ing  the  bubble  at  a  lower  angle  of  attack  (not  shown) 


begin  to  ‘open-up’  and  immediately  appear  above  the 
airfoil,  fringes  normal  to  the  surface  appear  near  the 
downstream  end  of  the  bubble.  This  can  be  seen  more 
clearly  at  q  =  15.5®.  By  a  =  15.75°,  these  fringes  ex¬ 
tend  to  about  1.5%  of  chord  above  the  airfoil  and 
nearly  to  0.2c  along  it.  The  inflexion  in  the  fringes 
very  close  to  the  airfoil  surface  suggests  that  there 
is  a  very  slight  reverse  flow.  In  contrast,  the  outer 
fringes  (outside  of  the  bubble)  proceed  only  in  one 
direction,  towards  the  trailing  edge.  The  region  of  re¬ 
verse  flow  is  less  than  0.005c  and  thus,  it  is  extremely 
difficult  to  detect  with  experimental  techniques  other 
than  surface  mounted  gages.  By  a  =  16.5°,  the  vor¬ 
tex  has  fully  developed  and  convected  to  0.2c.  This 
rapid  succession  of  events  in  a  very  small  angle  of  at¬ 
tack  range  is  typical  of  the  evolution  of  dynamic  stall 
at  compressible  Mach  numbers.  However,  the  inter¬ 
action  between  shocks  and  the  boundary  layer,  along 
with  the  events  of  bubble  bursting  and  vortex  for¬ 
mation  are  too  complicated  to  be  resolved  from  the 
images  processed  so  far. 

The  primary  issue  that  needs  to  be  determined 
is  the  origin  of  the  dynamic  stall  vortex.  Whether 
it  originates  independently  of  the  separation  bubble 
and  simply  pushes  the  back  end  of  the  bubble,  or  the 
increasing  angle  of  attack  makes  it  no  longer  possi¬ 
ble  for  the  bubble  to  remain  close  to  the  leading  edge 
and  moves  the  reattachment  point  rearward  allowing 
the  leading  edge  vorticity  to  coalesce,  is  yet  to  be 
found.  This  is  an  intriguing  issue,  especially  because 
the  leading  edge  flow  has  already  separated,  causing 
the  bubble,  at  a  very  low  angle  of  attack  (much  lower 
than  the  static  stall  angle).  A  plausible  description 
is  that  the  recirculation  region  in  the  bubble  becomes 
stronger  with  increase  of  angle  of  attack,  eventually 
forcing  a  breakdown  of  the  leading  edge  flow,  trig¬ 
gering  the  vortex  formation.  However,  a  much  more 
careful  analysis  is  needed  for  clarifying  the  issues. 

The  formation  of  a  separation  bubble  indicates 
that  the  airfoil  leading  edge  boundary  layer  transition 
plays  a  critical  role  in  dynamic  stall  occurrence.  If 
the  boundary  layer  is  modified  to  transition  before 
the  adverse  pressure  gradient  is  encountered,  it  is  to 
be  expected  that  dynamic  stall  and  the  associated 
events  occur  in  a  very  different  manner.  This  includes 
stall  onset,  the  details  of  the  vortex  formation  and  the 
duration  of  dynamic  lift  in  the  pitching  cycle.  This  is 
a  further  eispect  that  merits  a  more  detailed  study. 

F.  Trailing  Edge  Flow  Beyond  Deep  Stall 

The  flow  continues  to  show  interesting  features 
even  after  deep  stall  has  occurred  and  the  dynamic 
stall  vortex  is  shed.  In  Fig.  8,  the  complete  flow 
field  is  shown  for  just  one  experimental  flow  condi¬ 
tion  (M  =  0.45  and  o'*"  =  0.025).  For  this  case, 
the  dynamic  stall  angle  of  attack  (when  the  dynamic 
stall  vortex  leaves  the  airfoil)  is  17  degrees;  however, 
♦  he  nirfoil  continues  to  pl'ch  at  a  ,;onbtan*  rate  liii 
a  =  60°.  An  immediate  consequence  ot  vortex  shed¬ 
ding  under  these  conditions  is  an  increase  in  suction 
pressure  at  the  trailing  edge  on  the  lower  surface  of 
the  airfoil.  The  trailing  edge  suction  gradually  in¬ 
creases  during  the  dynamic  stall  process  as  the  airfoil 
pitches  from  Cp  Rs  0.0  at  q  =  0®  to  Cp  =  —0.504  at 
Q  =  17®.  However,  at  o  =  17.5®,  the  suction  pressure 
coefficient  jumps  to  -1.04.  Also,  a  counter  clockwise 
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vortex  begins  to  form  in  the  near-wake,  at  o  =  18“. 
The  birth/growth  of  this  vortex  puslies  the  separat¬ 
ing  shear  layer  from  the  leading  edge  outward.  At 
Q  =  20®,  this  vortex  is  shed  and  by  n  —  30°,  only 
the  separated  flow  can  be  seen  over  the  upper  sur¬ 
face.  The  same  number  of  fringes(26,  corresponding 
to  Cp  =  —1.08),  measured  from  the  stagnation  point 
to  this  point  on  the  upper  and  lower  surface  shear  lay¬ 
ers  seems  to  indicate  that  the  pressure  over  the  upper 
surface  is  nearly  constant.  However,  since  the  field 
of  view  is  limited  in  the  facility,  the  role  of  the  wake 
and  the  interaction  of  the  two  shear  layers  cannot  be 
determined  to  ascertain  this  fully.  As  the  pitch  up 
continues,  shear  layer  instabilities  develop,  which  roll 
up  into  vortices,  which  appear  at  a  =  ,50°  or  55°.  A 
large  trailing  edge  vortex  appears  from  the  lower  sur¬ 
face  at  a  =  50°,  which  rolls  up  towards  the  leading 
edge  shear  layer.  The  peak  suction  increases  slightly 
with  angle  of  attack,  even  though  the  airfoil  is  in  the 
deep  stall  state.  The  stagnation  point  moves  to  0.25 
chord  point  at  q  =  55°.  Similar  features  were  found 
for  other  Mach  numbers  as  well.  However,  this  pic¬ 
ture  is  significantly  different  from  the  low  Reynolds 
number  and  low  Mach  number  results  of  Walker  et 
al®,  who  observed  two  large  vortices  co-existing  on 
the  airfoil  suction  surface  as  dynamic  stall  progressed. 
This  led  to  a  much  different  airfoil  surface  pressure 
distribution  than  found  here.  For  all  the  cases  stud¬ 
ied  here,  the  trailing  edge  vortex  was  induced  as  the 
dynamic  stall  vortex  was  shed.  At  times  the  leading 
edge  separated  shear  layer  produced  a  large  vortex, 
resulting  a  double  vortex  pattern.  But,  it  then  had 
the  appearance  of  a  vortex  street. 

3.2.  Quantitative  Pressure  Field 

One  of  the  main  advantages  of  PDI  is  its  truly 
nonintrusive  way  of  yielding  the  pressure  field.  In 
a  rapidly  changing  flow  such  as  the  ramping  airfoil 
motion  under  consideration,  it  requires  several  real¬ 
izations  of  the  flow  -  i.e.  several  pitch  up  motions 
-  to  obtain  a  surface  pressure  distribution  by  stan¬ 
dard  techniques.  The  PDI  technique  not  only  pro¬ 
vides  the  surface  pressure  map,  but  also  the  global 
pressure  field  instantaneously  in  one  picture,  without 
any  history  effects  present  in  it.  Since  a  large  num¬ 
ber  of  interferograms  were  obtained  at  close  intervals 
in  angle  of  attack,  a  finer  set  of  instantaneous  pres¬ 
sure  data  is  in  hand  now.  The  data  to  be  reported 
are  the  first  global  mapping  pressure  coefficients  that 
have  been  obtained  for  a  transiently  pitching  airfoil 
under  compressibility  conditions  using  PDI. 

A.  Effect  of  Mach  Number  on  Peak  Suction  Pressure 

Fig.  9  shows  the  peak  suction  pressure  coefficient 
obtained  over  the  airfoil  at  different  Mach  numbers 
for  different  pitch  rates.  The  absolute  value  of  Cp 
in  steady  flow  is  lower  than  in  unsteady  flows  for  all 
Mach  numbers  tested.  Of  significance  is  the  relative 
independence  of  the  peak  suction  pressure  coefficient 
of  the  pitch  rate  for  all  Mach  numbers.  Also,  the  peak 
suction  decreases  dramatically  with  increasing  Mach 
number.  For  the  lowest  Mach  number  of  0.2,  the  value 
of  Cp  =  —7.1  corresponds  to  a  local  Mach  number  of 
0.58,  indicating  that  compressibility  conditions  play 


a  role  in  the  process.  At  M  =  0.45.  the  sonic  limit 
is  Cp  =  —2.78  and  with  the  suction  peak  at  Cp  = 
—3.6,  the  local  flow  is  supersonic  and  thus  is  strongly 
compressible.  It  is  possible  that  the  formation  of  the 
bubble  and  the  strong  local  compressible  conditions 
have  both  affected  the  viscous/inviscid  interactions  in 
the  flow  and  thus,  limited  the  suction  pressure  from 
increasing  as  the  pitch  rate  is  increased. 

B.  Effect  of  Pitch  Rate 

In  Fig.  10,  the  development  of  the  leading 
edge  pressure  distribution  over  the  first  5%  of  air¬ 
foil  chord  is  compared  for  M  =  0.3  for  steady  flow; 
Q+  =  0.035  and  a"*"  =  0.04  at  a  =  10  degrees.  It 
is  clear  that  steady  flow  develops  a  stronger  suction 
peak  than  the  unsteady  flow  Ccises,  with  a  maximum 
Cp  of  -3.6;  the  suction  peak  is  located  at  x/c  =  0.012. 
Tne  adverse  pressure  gradient  is  slower  to  develop  at 
the  higher  pitch  rates;  in  fact,  the  suction  level  has 
reached  only  a  value  of  Cp  =  —2.95  at  10  degrees 
angle  of  attack.  The  pressure  distribution  for  the  un¬ 
steady  cases  is  flatter  than  in  the  steady  flow  case 
as  observed  from  the  interferograms.  Even  though 
all  three  cases  show  a  separation  bubble,  as  observed 
on  the  interferograms,  tb^  effect  of  increasing  pitch 
rate  is  clearly  to  delay  the  flow  development  over  the 
leading  edge,  which  consequently  leads  to  delay  in  dy¬ 
namic  stall  occurrence. 

Fig.  11  shows  the  maximum  suction  pressure  co¬ 
efficient  at  different  angles  of  attack  for  M  =  0.4.  In 
it,  steady  flow,  and  unsteady  flows  at  =  0.02  and 
a"*"  =  0.035  are  compared.  The  steady  flow  stalls  at 
a  =  10.8°;  at  higher  angles  the  leading  edge  suction 
falls  rapidly.  At  o+  =  0.02  the  effect  of  decreased 
adverse  pressure  gradient  as  seen  in  Fig.  8  is  seen  at 
all  angles  of  attack;  the  suction  levels  for  o'*'  =  0.035 
lag  those  of  even  =  0.02.  The  rounding  of  the 
curves  at  the  top  corresponds  to  the  situation  when 
the  dynamic  stall  vortex  is  in  the  formative  stages. 
During  this  stage  (q  =  12°  —  13°  for  a"*"  =  0.02 
and  a  =  13°  —  13.5°  for  o+  =  0.035)  the  airfoil 
suction  is  the  maximum.  Once  the  vortex  begins  to 
convect,  the  leading  edge  flow  slows  down  and  the  suc¬ 
tion  is  steadily  lost,  even  though  the  angle  of  attack 
is  still  increasing.  Even  though  there  is  experimental 
scatter,  it  is  also  clear  that  dynamically  pitching  air¬ 
foils  can  withstand  larger  suction  peaks  than  steady 
airfoils  prior  to  stall  and  withstand  higher  flow  gra¬ 
dients,  before  stall  occurs.  It  is  also  of  value  to  note 
that  the  loss  of  suction  peak  in  the  dynamic  case  does 
not  mean  loss  of  lift,  but  only  the  initiation  of  the  dy¬ 
namic  stall  vortex. 

Fig.  12  shows  the  pressure  distribution  at  M  = 
0.3  at  Q  =  10°  in  steady  flow  and  q  =  11°  in  unsteady 
flow  at  a"*"  =  0.035.  The  nearly  identical  pressure 
distributions  suggest  that  unsteadiness  essentially  has 
introduced  a  one  degree  lag  in  the  effective  angle  of  at¬ 
tack,  including  the  formation  of  the  separation  bubble 
(see  next  section),  supporting  the  conclusions  drawn 
from  Fig.  10  and  Fig.  11.  As  much  as  two  degrees  of 
decrease  in  the  effective  angle  of  attack  was  observed 
for  oscillating  airfoils  by  Carr  et  al’^.  However,  as  can 
be  expected,  the  amount  of  this  benefit  is  determined 
by  the  experimental  conditions,  especially  the  pitch 
rate  and  Mach  number. 
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C.  Airfoil  Pressure  Distributions 

The  airfoil  pressure  distributions  presented  in 
Fig.  13  were  obtained  for  different  Mach  numbers 
and  pitch  rates  using  the  method  described  earlier. 
Only  one  typical  case  is  discussed  below  for  M  =  0.3, 
a"*"  =  0.035  at  different  angles  of  attack  as  seen  in 
Fig.  13.  For  this  case,  Fig.  13a  shows  that  at  a  =  5", 
the  pressure  distribution  is  smooth,  with  the  suction 
peak  of  Cp  =  —1.60  at  x/c  =  0.025.  As  the  airfoil 
pitches  to  7  degrees,  Cp  increases  to  -2.21.  This  is  ac¬ 
companied  by  the  upstream  movement  of  the  suction 
peak  to  x/c  =  0.015.  Along  with  this,  the  stagna¬ 
tion  point  moves  from  x/c  =  0.01  to  x/c  =  0.05.  As 
the  airfoil  pitches  to  higher  angles  of  attack,  further 
movement  of  the  stagnation  point  is  difficult  to  de¬ 
tect,  until  larger  angles  are  reached.  At  o  =  8®,  a 
separation  bubble  forms.  The  pressure  plateau  fol¬ 
lowing  the  steep  drop  in  the  pressure  distribution  im¬ 
mediately  after  the  suction  peak  at  a  =  9°  confirms 
this.  The  bubble  extends  to  x/c  «  0.05.  The  suction 
level  and  the  adverse  pressure  gradient  both  continue 
to  increase  as  the  airfoil  angle  of  attack  is  increased. 
However,  Fig.  13b  shows  that  the  pressure  distribu¬ 
tions  at  a  =  13°  and  a  =  14°  are  only  incrementally 
different  from  that  at  a  =  12°,  indicating  that  the 
lift  and  thus,  the  circulation  over  the  airfoil  have  not 
changed  significantly  in  this  range  of  pitching.  It  is 
also  the  range  when  the  dynamic  stall  vortex  is  initi¬ 
ated  as  determined  by  visual  analysis  of  the  interfer- 
ogram  images.  The  suction  peak,  however,  increases 
with  angle  of  attack  and  reaches  its  maximum  value  of 
-5.8  when  the  first  imprint  of  dynamic  stall  is  clearly 
discernible  in  the  images.  At  this  instant,  it  appears 
that  the  vortex  becomes  fully  organized  and  then  it 
begins  to  convect.  The  movement  of  the  vortex  re¬ 
sults  in  a  slight  drop  of  airfoil  suction.  However,  the 
pressure  distribution  changes  to  the  nearly  flat-top 
shape  which  extends  to  about  x/c  =  0.3.  This  in¬ 
creasing  region  of  low  pressure  over  the  upper  surface 
is  the  cause  of  the  vortical  lift  seen  in  the  flow.  Fur¬ 
ther  pitching  of  the  airfoil  continues  to  decrease  the 
suction  levels,  while  spreading  its  extent  more  over 
the  airfoil,  until  it  is  shed  past  the  trailing  edge. 

4.  Concluding  Remarks 

A  nonintrusive  study  of  the  compressibility  effects 
on  dynamic  stall  of  a  transiently  pitching  airfoil  has 
been  conducted.  The  study  shows  the  following: 

1.  Dynamic  stall  occurs  rapidly  as  the  separa¬ 
tion  bubble  bursts.  This  points  to  the  need  for  bet¬ 
ter  experimental  and  computational  modelling  of  the 
boundary  layer  transition  process. 

2.  Multiple  shocks  form  over  the  airfoil  for  a  free 
stream  Mach  number  of  0.4  and  higher.  The  role  of 
the  shocks  in  altering  the  dynamic  stall  process  needs 
a  careful  study. 

3.  Compressibility  lowers  the  angle  of  attack  at 
which  stall  occurs,  primarily  due  to  unfavorable  in¬ 
teraction  with  the  airfoil  boundary  layer. 

4.  The  duration  of  dynamic  lift  lasts  over  a  con¬ 
siderably  smaller  angle  of  attack  range  when  com¬ 
pressibility  effects  are  strong. 

5.  There  is  marked  differences  in  the  vortex 


structure  between  the  compressible  and  incompress¬ 
ible  cases. 

6.  Increcising  Mach  number  decreases  the  max¬ 
imum  Cp  developed  by  the  airfoil.  However,  for  a 
given  Mach  number  it  is  relatively  independent  of  the 
pitch  rate. 

7.  Interferogram  images  confirm  the  low  pressure 
region  over  the  airfoil  upper  surface  and  support  the 
earlier  result  that  dynamic  lift  is  generated  for  com¬ 
pressible  conditions. 

8.  Dynamic  motion  causes  a  lag  in  the  pressure 
field  development  contributing  to  stall  delay. 

Further  analysis  of  the  leading  edge  flow  images 
is  underway. 
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Fig.  7.  Details  of  Dy 
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Fig.  10.  Effect  of  Pitch  Rate  on  Pres.sure  Development  over  a  Transiently  Pitching  Airfoil,  M  =  O  .l 
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Fig.  11.  Effect  of  Pitch  Rate  on  Peak  Suction  Pressure  Coefficient  ov.;r  a  Transiently  Pitching 

Airfoil,  M  =  0.4. 
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Fig.  12.  Comparison  of  Steady  and  Unsteady  Flows,  M  =  0.3. 
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